Coherent control of atomic spin currents in a double well 
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We propose a method for controlling the atomic currents of a two-component Bose-Einstein con- 
densate in a double well by applying an external field to the atoms in one of the potential wells. We 
study the ground-state properties of the system and show that the directions of spin currents and 
net-particle tunneling can be manipulated by adiabatically varying the coupling strength between 
the atoms and the field. This system can be used for studying spin and tunneling phenomena across 
a wide range of interaction parameters. In addition, spin-squeezed states can be generated. It is 
useful for quantum information processing and quantum metrology. 
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I. INTRODUCTION 

Spin and tunneling phenomena are of fundamental in- 
terests in understanding quantum behaviors of particles. 
They are also important in the applications using, solid- 
state devices such as sensors and data storage [Ij. In 
addition, manipulation of the quantum state of a single 
spin is essential in implementing quantum information 
processing Q. 

Recently, the tunneling dynamics of ultracold atoms 
has been observed in a double well Q and optical lattices 
[H], where the experimental parameters can be widely 
tuned. Moreover, high-fidelity single-spin detection of 
an atom has been realized in an optical lattice [f| 
and atom-chip Q, respectively. Such sophisticated tech- 
niques of manipulating ultracold atoms open up the pos- 
sibilities for the study of intriguing quantum phenomena 
never possible before. 

In this paper, we propose a method to control the tun- 
neling dynamics of a two-component Bose-Einstein con- 
densate (BEC) in a double well H by applying an exter- 
nal field to the atoms in one of the potential wells. In 
fact, the methods for controlling of tunneling in a double- 
well BEC have been suggested by driving the double- well 
potentials Q and by applying a symmetry-breaking field 
10], respectively. Besides, a number of methods for ma- 
nipulating the atomic motions in an an optical lattice 
have been proposed such as using external fields for vi- 
brational transitions between adjacent sites [ll|, tilting 
the lattice potential [jjl and periodic modulation of the 
lattice parameters (l3j | . 

Here we show that the spin and tunneling dynamics of 
atoms can be manipulated by adiabatically changing the 
coupling strength of the field. This approach can be used 
for studying the spin and tunneling related phenomena 
in a controllable manner. For example, the directions of 
"spin currents" such as parallel- and counter-flows can 
be controlled by appropriately adjusting the interaction 
parameters. Here the "spin currents" refer to the two 
atomic currents of a condensate of 87 Rb atoms with the 
two different hyperfme levels [lij ]. Apart from changing 



the internal states of atoms, the external field gives rise 
to net-particle tunneling. It is different to the situation of 
co-tunneling of the two component condensates in a dou- 
ble well 0, llo| in which the number difference of atoms 
between the wells is equal during the tunneling process. 

In addition, the tunnel behaviors are totally different in 
the limits of weak and strong atomic interactions. In the 
regime of weak atomic interactions, the atoms smoothly 
tunnel through the other well. On the contrary, the dis- 
crete steps of tunneling are shown in the limit of strong 
atomic interactions. The coherent control of single-atom 
tunneling can thus be achieved. This method can be uti- 
lized for atomic transport. 

Furthermore, we investigate the production of spin 
squeezing [lfl [I?) in this system. The occurrence of spin- 
squeezing can indicate multi-particle entanglement [18]. 
We show that spin-squeezed states can be dynamically 
generated by slowly changing the coupling strength of 
the field. This can be used for preparing entangled states 
and quantum metrology fljjj . 

This paper is organized as follows: In Sec. II, we in- 
troduce the system of a two-component BEC in a double 
well, where the atoms in the left potential well are cou- 
pled to the laser field. In Sec. Ill, we study the ground- 
state properties of the system. In Sec. IV, we propose a 
method to adiabatically transfer the atoms to the other 
well. In Sec. V, we investigate the generation of multi- 
particle entanglement using this method. We provide a 
summary in Sec. VI. In the Appendix, we derive an ef- 
fective Hamiltonian describing the coupling between the 
two internal states with the lasers. 



II. SYSTEM 

We consider a condensate of 87 Rb atoms with two hy- 
perfme levels \e) = \F = 2,mp = 1) and Jo) = \F = 
l,mp = —1) of the hSi/2 ground state [1J| confined 
in a symmetric double- well potential An external 
field is applied to the atoms in one of the potential 
wells. The schematic of the system is shown in Fig.[TJa). 
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FIG. 1. (Color online) (a) Schematic of a two-component 
Bose-Einstein condensate in a double well. An external field 
is applied to the atoms in the left potential well, (b) Energy 
levels for the atoms. The two internal states |e) and \g) are 
coupled via the upper state \r). 

This system can be described by the total Hamiltonian 
i/tot = Hq + H\ , where Hq and H i are the Hamiltonians 
describing the external and internal degrees of atoms. 

We adopt the two-mode approximation to describe the 
atoms in deep potential wells Q. Since the scatter- 
ing lengths of the different hyperfine states of 87 Rb are 
very similar [2(|, we assume that the intra- and inter- 
component interactions are nearly the same. The Hamil- 
tonian Hq can be written as (ft = 1) 

H = -J{e f L e R + e R e L + g\g R + g R g L ) 

+U[{n eL + n gL f + (n eR + n gR f], (1) 

where (Xl(cxr) and n aL (n aR ) are annihilation operator 
of an atom and number operator in the left (right) po- 
tential well, for a — e, g. The total number N of atoms 
is conserved in this system. To ensure the validity of 
two-mode approximation, we assume that the trapping 
energy is much larger than the atomic interaction energy 
SUll, i-e., wo 3> UN, where ujq is the effective trapping 
frequency of the potential well. 

We can make a rough estimation of the experimen- 
tal parameters within the two-mode approximation. For 
87 Rb, the scattering length a is about 5 nm. We take the 
frequencies u y and uj z of the transervse trapping poten- 
tials to be about 2irxl kHz. Indeed, the barrier height 
and the separation between the two potential minima 
can be varied [22J by appropriately splitting the poten- 
tial. We consider that the barrier height Vb can be tuned 
22] from 27r x 50 to 27r x 250 Hz and the separation 2xq 
between two w ells are 2 a nd 4 /an. The effective fre- 
quency cj = \/SVb/ { mx o) i s related to the barrier 
height Vb and the separation 2xo. This gives the effective 
frequency w ranging from 2-7T x 220 to 2ir x 350 Hz, and 
the interaction strength U ranges between 2ir x 2.7 and 
27rx3.4 Hz. From these estimations, the number N of 
atoms must be less than 100 to maintain the validity of 
two- mode approximation [ij |2lj. We also estimate the 
ratio U I J of the atomic strength and tunneling strength 
which ranges from 0.05 to 300. 

Without loss of generality, we consider an external field 
to be applied to the atoms in the left potential well. The 
two internal states can be coupled via the upper transi- 
tion of the D2 line of 87 Rb 23j by using the two laser 



beams with different circular polarizations as shown in 
Fig- Hlb). This upper state \r) can be adiabatically elim- 
inated due to the large detuning. In the Appendix A, we 
derive an effective Hamiltonian for describing the inter- 
action between the two internal states and the lasers. In 
the interaction picture, the Hamiltonian Hi is given by 




H\ = A(n eL + n eR )+m(elg L +n.c), (2) 

where A is the detuning between the atomic transition 
(\r) and |e)) and the laser field, and il is the effective 
coupling strength between the atoms and external field. 
A tightly focused laser can be applied to the atoms in one 
of the potential well. In fact, a tightly focused laser has 
been used for addressing a single atom in an optical lat- 
tice @, where the full- width at half- maximum (FWHM) 
of the diameter of laser beam is within 1 fim [6\ . The sep- 
aration between the two potential wells is about several 
/im in the experiment Therefore, the effects of the 
external lasers on the atoms in the other well are small. 



III. GROUND-STATE PROPERTIES OF THE 
COUPLED ATOM-LASER SYSTEM 

Now we study how the ground state properties of the 
BEC affected by the local external field. In Fig. [H we 
plot the population differences {{n a i) — (n aR )) versus the 
coupling strengths f2 for the atoms in the two different 
internal states \a) and a — e,g. For the cases of even 
number N of atoms, there is an equal number of atoms 
in the two wells in the absence of the external field, i.e., 
£7 = 0. The external field causes the energy bias between 
the two wells. Thus, the population difference becomes 
larger when the coupling strength increases. 

Moreover, the system exhibits totally different behav- 
iors in the regimes of weak and strong atomic inter- 
actions. For weak atomic interactions, the population 
differences smoothly vary as a function of the coupling 
strength as shown in Figs. HJa) and (c). Also, a larger 
number of atoms are in the state \g) due to the larger 
detuning A. 

In Figs.[UJb) and (d), we plot the population differences 
versus the coupling strengths fi in the regime of strong 
atomic interactions. We can see that discrete steps of 
the population difference for atoms in the state \g) are 
shown when the coupling strength increases. However, 
the discrete feature is not obvious for the atoms in the 
state |e). In Fig. [UJb), the atoms, in the two different 
internal states, distribute in the opposite potential wells 
for small ft. When ft becomes larger, both component 
of atoms populate in the left potential well. This result 
shows that the population difference of atoms in the two 
internal states depends on the coupling strength and also 
the detuning between the atoms and field. 

To proceed, we investigate the relationship between 
the total population difference of atoms ((n e £ + n g L — 
n e R ~ n gR )) and the coupling strength Q. In Fig. [31 we 
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FIG. 2. (Color online) Population differences are plotted 
as a function of Sty J for the ground state of the system and 
N = 8. The parameters are shown: (a) and (b) A = 20J; (c) 
and (d) A = 200J. The ground state \g) and excited state 
|e) of atoms are denoted by blue-solid and red-dotted lines, 
respectively. 



plot the total population differences as a function of the 
coupling strength SI for different strengths U of atomic 
interactions. The external field leads to the population 
imbalance between two wells in both regimes of weak 
and strong atomic interactions. For weak atomic interac- 
tions, the population difference smoothly increases with 
SI. When the atomic interactions become strong, the dis- 
crete steps of population differences are shown in Fig. [3] 
The sharper discrete steps are shown for larger U. In this 
case, a single atom is only allowed to tunnel through the 
other well for the specific coupling strengths. 



A. Tunneling condition in the limit of strong 
atomic interactions 

We now discuss the tunneling condition in the limit of 
strong atomic interactions. Since the tunnel couplings 
are negligible in this regime, the numbers of atoms in 
the two wells are conserved. We assume that there are 
N/2 + n and N/2 — n atoms in the left and right wells. 
For convenience, we define the angular momentum oper- 
ators: S jx = (gjeJj + ejgt)/2, S jy = (g s e^ - e j g} j )/2i and 
Sj Z = (eje 3 — g\gj)/2, and j — L,R. To diagonalizc the 



Hamiltonian H to t, we apply the transformation as 



Slx = cos8S' Lx -sm9S' l 
S Lz = cos9S' Lz + sm.eSi 



J Lz"> 



Lx ' 



(3) 
(4) 



By setting the term, A sin 8 + 2S1 cos 9, to zero, the total 
Hamiltonian can be diagonalized as 



H'tot 



AS Rz + hy/A 2 + 4fi 2 S£ 2 + U(N 2 /2 + 2n 2 ) 
+HAN/2, 
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FIG. 3. (Color online) Plot of the total population difference 
versus SI/ J for the ground state with the different strengths 
U of atomic interactions. The parameters are used: N = 8 
and A = 20J. The values of S7„ in Eq. ([6]) are marked with 
the red empty circles. 



and its ground-state energy is given by 



E, 



-A(JV/2-n)/2- 
+2Un 2 + UN 2 /2- 



AQ 2 (N/2 + n)/2 



AN/2. 



(6) 



The tunneling of atoms occurs when the energies E~* 
and are equal to each other. In this situation, a 

single atom tunnels through the other well. It is anal- 
ogous to the resonant tunneling in quantum dots due 
to the Coulomb blockade 0. By setting E^ = E^_ x in 
Eq. ([6]) , the tunneling condition for the coupling strength 
fi n can be obtained as 



= \{W{2n-l) 



A] 2 - A 2 p. 



(7) 



In Fig. [31 the values of f2„ are marked with the red 
empty circles. This shows that the tunneling condition 
in Eq. Q agrees with the exact numerical solution. 



IV. ADIABATIC TRANSPORT 

We have studied the ground-state properties of the 
coupled atom-laser system in the previous section. We 
have shown that the directions of atomic currents and the 
population difference depend on the coupling strength 
n. Now we study the adiabatic transport of atoms by 
slowly increasing the coupling strength tt. According to 
the adiabatic theorem [2(|, the system can evolve as its 
instantaneous ground state if the changing rate of the 
parameter Q is sufficiently slow. Therefore, this method 
can be used for adiabatically transferring the atoms to 
the other potential well. 

We consider the coupling strength Q(t) as a linear func- 
tion of time t, i.e., 



tt(t) = vt, 



(8) 



(5) 



where v is a positive number. Here the detuning A is kept 
as a constant during the time evolution. Initially, the 
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FIG. 4. (Color online) The population differences for the 
two different component condensates are plotted against the 
time Jt for U — 0.1J and 10J in (a) and (b). The atoms in 
ground and excited states are denoted by blue-solid and red- 
dashed lines, respectively, (c) Plot of the total population 
difference versus the time Jt for different strengths of atomic 
interactions. The parameters are used: N = 8, A = 20J and 
v = 2.5J. 



FIG. 5. (Color online) The total population differences are 
plotted versus the time Jt for the different detunings A = 20 J 
and 200 J in (a) and (b), respectively. The different rates of 
change are denoted: v = 10 J (black-solid line), v = 20 J 
(blue-dashed line), v — 30 J (red-dotted line) and v = 40 J 
(green-dashed-dotted line), respectively. The parameters are 
used: N = 8 and U = 10 J. 



A. Efficiency of the population difference 



system is prepared in its ground state of the Hamiltonian 
Ho in Eq. (JTJ and setting O = 0. The coupling strength 
f2(t) in Eq. ((5J) is adiabatically increased. In Figs. Ufa) 
and (b) , the population differences are plotted versus the 
time for the atoms in the different internal states. We 
can see that the two results in Figs. [5] and 0] reach a good 
agreement. This shows that the tunneling dynamics of 
atoms can be controlled by using an external field. 

In the regime of weak atomic interactions, the two dif- 
ferent component condensates smoothly tunnel through 
the other well in the same direction in Fig. 0|a). In the 
opposite interaction limit, the discrete steps are shown in 
Fig. Bib). Note that the counter-flow is shown in a short 
time in figure (b) . The flows of atomic spin currents be- 
come parallel afterward. This shows that the direction of 
spin flows can be controlled by appropriately adjusting 
the parameters SI and A. 

Then, we study the time evolution of total population 
difference of atoms. In Fig. EJc), we examine the tun- 
neling dynamics for a wide range of interaction param- 
eters. In both regimes, the atoms tunnel to the left po- 
tential well when the coupling strength slowly increases 
with time. The total population differences smoothly 
increases as a function of time for weak atomic interac- 
tions. The discrete steps of the tunneling are shown when 
the strength of atomic interactions becomes strong. The 
single-atom tunneling can thus be achieved in the limit 
of strong atomic interactions. 



Next, we investigate the efficiency of the population 
transfer by increasing the coupling strength of the exter- 
nal field. In Figs.OJa) and (b), the population differences 
are plotted versus the time for the two different detun- 
ings A = 20 J and 200 J. The full population transfer can 
be achieved if the parameters v are small enough in both 
cases. When the parameters v become larger, the rates of 
population transfer increase. However, the smaller pop- 
ulation of atoms can be transferred in both cases. It is 
because they have gone beyond the adiabatic limit. By 
comparing the figuresjnia) and (b), we find that the larger 
numbers of atoms can be transported with the same rate 
v of change for the case using a larger detuning. 

B. Effects of the coupling between the laser and 
the atoms in the neighboring well 

Since the separation between the two wells is small 
[22} , the lasers may also couple the atoms in the other 
potential well. Now we examine a very small coupling 
between the laser and the atoms in the right potential 
well. The Hamiltonian, describes the coupling between 
the laser and the atoms in the right potential well, can 
be written as 

H? = n'(e R g R + g R e R ), (9) 

where ST is the coupling strength between the laser and 
atoms. Here we consider that the coupling strength SI' 
linearly increases with time as 

n'(t) = v't. (10) 
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FIG. 6. (Color online) Time evolution of the total population 
difference. The effect of a small coupling between the laser 
and the atoms in the right potential well is investigated. In 
(a), the total population differences are plotted versus the 
time Jt for different rates of changing v' and N = 8. In (b), 
the total population differences are plotted versus the time Jt 
for different numbers of atoms and v' — J. The parameters 
are used: U = 10 J, A = 20J and v = 10J. 



We assume that the parameter v' is much smaller than v 
in Eq. ©. 

In Fig. [HJa), the total population differences are plot- 
ted versus the time for the different rates of changing v' . 
The atoms can be efficiently transferred from the left to 
right potential well. The results show slightly different 
tunneling behaviors for small v' . We then examine the 
population transfer with small v' for larger number of 
atoms. In Fig. Efb), we plot the total population dif- 
ference as a function of time for different numbers of 
atoms. The discrete steps of tunneling can be clearly 
shown. This shows that this method works even if there 
is a very small coupling between the laser and the atoms 
in the neighboring well. 



V. MULTI-PARTICLE ENTANGLEMENT 

Having discussed the tunneling dynamics of atoms, we 
study the generation of spin-squeezed states by adiabat- 
ically changing the coupling strength of the field. To 
indicate the occurrence of spin squeezing, a parameter 
£ 2 can be defined as [llj 



N(AS ni f 



{S n 



:) 2 



r 



(ii) 



where rii is the i-th component of an angular momentum 
system, and i = 1,2 and 3. If £ is less than one, then 
the system is said to be spin-squeezed |17| . In addition, 
the parameter £ can be used for indicating multi-particle 
entanglement [la ]. 

Let us define the angular momentum operators as: 
S x = (4,e fi + g\g R + e R e L + g R g L )/2, S y = (e{e R + 



Zs0.5 




FIG. 7. (Color online) The spin-squeezing parameters £ 2 are 
plotted versus the time Jt for different number N of atoms. 
The parameters are used: A = 200J, U = 10J and v = 10J. 



-9l9l- 



-e R e R - 



e R e L-9 R 9L)/2i, and S z = (e\e R - 
g R g R )/2. The angular momentum operators obey the 
standard commutation rule. We study the parameter as 



£ 2 = 



N(AS Z ) 2 
(S x ) 2 



(12) 



Physically speaking, the quantity (AS Z ) 2 is the variance 
of the total population difference of atoms between the 
wells, and (S x ) is the sum of the phase coherences be- 
tween the two wells for the two component condensates. 

In Fig. [71 we plot the spin-squeezing parameter £ 2 ver- 
sus the time for different number of atoms. Initially, the 
parameter £ 2 is below one when fl — 0. This means that 
the initial ground state is a spin-squeezed state. As 51 
in Eq. ([7J slowly increases with time, spin squeezing can 
be dynamically produced. This result indicates that the 
system is spin-squeezed for a wide range of f2. Besides, 
a higher degree of spin squeezing can be produced with 
a larger number of atoms N. 



VI. CONCLUSIONS 

In summary, we have studied how the ground state 
of a two-component condensate in a double well affected 
by a local external field. We have shown that the flows 
of spin currents and particle-tunneling dynamics can be 
controlled by slowly varying the coupling strength of the 
external field and appropriately adjusting the detuning. 
This can be used for studying spin and tunneling phe- 
nomena and also the potential applications of atomic 
devices in atomtronics [27| . In addition, spin-squeezed 
states can be dynamical generated. It is an important 
resource for precision measurement 
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Appendix A: Derivation of the effective Hamiltonian 
for coupling between two internal states and the 
lasers 



The two internal states of an atom can be coupled via 
pumping to the upper state by two lasers as shown in 
Fig. [Tfb) . The Hamiltonian, describes the interaction 
between the two internal states and the lasers, can be 
written as 

H = uj rg \r)(r\ +u; e9 |e)(e| + fi s [exp (-iu g t)\r)(g\ + H.c] 
+^ e [exp (-iu) e t)\r)(e\ + H.c], (Al) 

By performing the unitary transformation as, 

U{t) = cxp{-i[Lo g \r)(r\ + {oj g + w e )|e)<e|]t}, (A2) 

the Hamiltonian can be transformed as [24{ 

H 1 = iU^U + HU, (A3) 
= A r \r){r\ + A e \e){e\ + {n g \r)(g\ + n e \r)(e\ + H.c), 

(A4) 



where A r 



iO rq 



uj g and A e 



We write the state |^(t)) as, 

\*(t)) = c g \g)+c e \e)+c r \r). (A5) 
From the Schrodinger equation, we can obtain 



iCg QgC ri 



= A T 



(A6) 
(A7) 
(A8) 



We assume that £l e and £l g are real numbers. 

If the detuning A is much greater than the coupling 
strengths tt e ,£l g , then the upper state can be adiabati- 
cally eliminated. We can obtain 



tig 

a/ 9 



A r 



(A9) 



The effective Hamiltonian can thus be obtained as 



H cS = -^\g)(g\ + (A e - ^) |e)(e| + n(\g)(e\ +H.c), 

(A10) 

where fl = —il e fl g /A r is the effective coupling strength 
between the two internal states \g) and |e). We assume 
that fig and £l e are approximately equal and let A e = A 
in Eq. ©. 
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